Belle Preprint 2008-17 
KEK Preprint 2008-12 



O 

(N 

C 
3 



Time-dependent CP Asymmetries in B° — > Kgp ^ Decays 

J. Li, 6 I. Adachi, 7 K. Arinstcin, 1 T. Aushev, 17 * 12 A. M. Bakich, 37 V. Balagura, 12 I. Bedny, 1 V. Bhardwaj, 32 
U. Bitenc, 13 A. Bozek, 26 M. Bracko, 19 * 13 T. E. Browder, 6 P. Chang, 25 Y. Chao, 25 A. Chen, 23 B. G. Cheon, 5 
R. Chistov, 12 Y. Choi, 36 J. Dalseno, 7 A. Drutskoy, 3 S. Eidelman, 1 N. Gabyshev, 1 H. Ha, 15 K. Hara, 21 
Y. Hasegawa, 35 H. Hayashii, 22 M. Hazumi, 7 D. Heffernan, 31 Y. Hoshi, 40 W.-S. Hou, 25 H. J. Hyun, 16 T. Iijima, 21 
A. Ishikawa, 33 R. Itoh, 7 M. Iwasaki, 42 Y. Iwasaki, 7 N. J. Joshi, 38 D. H. Kah, 16 J. H. Kang, 46 H. Kawai, 2 
T. Kawasaki, 28 H. Kichimi, 7 Y. I. Kim, 16 Y. J. Kim, 4 K. Kinoshita, 3 S. Korpar, 19 ' 13 P. Krizan, 18 - 13 P. Krokovny, 7 
R. Kumar, 32 A. Kuzmin, 1 S.-H. Kycong, 46 C. Liu, 34 Y. Liu 4 A. Matyja, 26 S. McOnie, 37 T. Medvedeva, 12 
K. Miyabayashi, 22 H. Miyake, 31 H. Miyata, 28 G. R. Moloney, 20 Y. Nagasaka, 8 M. Nakao, 7 Z. Natkaniec, 26 
S. Nishida, 7 O. Nitoh, 44 S. Ogawa, 39 T. Ohshima, 21 S. Okuno, 14 S. L. Olsen, 6 ' 9 H. Ozaki, 7 G. Pakhlova, 12 

C. W. Park, 36 H. Park, 16 H. K. Park, 16 K. S. Park, 36 L. S. Peak, 37 R. Pestotnik, 13 L. E. Piiloncn, 45 
H. Sahoo, 6 Y. Sakai, 7 O. Schneider, 17 C. Schwanda, 10 A. Sekiya, 22 K. Senyo, 21 M. Shapkin, 11 J.-G. Shiu, 25 
B. Shwartz, 1 A. Sokolov, 11 A. Somov, 3 S. Stanic, 29 M. Staric, 13 T. Sumiyoshi, 43 M. Tanaka, 7 G. N. Taylor, 20 
Y. Teramoto, 30 I. Tikhomirov, 12 K. Trabelsi, 7 T. Tsuboyama, 7 S. Uehara, 7 T. Uglov, 12 Y. Unno, 5 S. Uno, 7 
Os| i P. Urquijo, 20 Y. Ushiroda, 7 Y. Usov, 1 G. Varner, 6 K. E. Varvell, 37 K. Vervink, 17 A. Vinokurova, 1 

C. H. Wang, 24 P. Wang, 9 X. L. Wang, 9 Y. Watanabe, 14 E. Won, 15 B. D. Yabsley, 37 H. Yamamoto, 41 
Y. Yamashita, 27 M. Yamauchi, 7 Z. P. Zhang, 34 V. Zhilich, 1 T. Zivko, 13 A. Zupanc, 13 and O. Zyukova 1 

(The Belle Collaboration) 

Budker Institute of Nuclear Physics, Novosibirsk 
Chiba University, Chiba 
^ . 3 University of Cincinnati, Cincinnati, Ohio 45221 

J * The Graduate University for Advanced Studies, Hayama 
5 Hanyang University, Seoul 
6 University of Hawaii, Honolulu, Hawaii 96822 
^ ' 7 High Energy Accelerator Research Organization (KEK), Tsukuba 

, 8 Hiroshima Institute of Technology, Hiroshima 

00 ■ 9 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 

Institute of High Energy Physics, Vienna 
11 Institute of High Energy Physics, Protvino 
\Q • 12 Institute for Theoretical and Experimental Physics, Moscow 

' 13 J. Stefan Institute, Ljubljana 

OO , 1J> Kanagawa University, Yokohama 

. 15 Korea University, Seoul 

lb Kyungpook National University, Taegu 
11 E cole Polytechnique Eederale de Lausanne (EPFL), Lausanne 
■ ls Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana 

' 19 University of Maribor, Maribor 

University of Melbourne, School of Physics, Victoria 3010 
21 Nagoya University, Nagoya 
22 Nara Women's University, Nara 
2:1 National Central University, Chung-li 
2Jl National United University, Miao Li 
25 Department of Physics, National Taiwan University, Taipei 
26 H. Niewodniczanski Institute of Nuclear Physics, Krakow 
21 Nippon Dental University, Niigata 
28 Niigata University, Niigata 
29 University of Nova Gorica, Nova Gorica 
30 Osaka City University, Osaka 

31 Osaka University, Osaka 
32 Panjab University, Chandigarh 
33 Saga University, Saga 
3J * University of Science and Technology of China, Hefei 
35 Shinshu University, Nagano 
36 Sungkyunkwan University, Suwon 



37 University of Sydney, Sydney, New South Wales 



Typeset by REVTeX 



2 



38 Tata Institute of Fundamental Research, Mumbai 
39 Toho University, Funabashi 
J * u Tohoku Gakuin University, Tagajo 
Jfl Tohoku University, Sendai 
Jf2 Department of Physics, University of Tokyo, Tokyo 

JfS Tokyo Metropolitan University, Tokyo 
Tokyo University of Agriculture and Technology, Tokyo 
45 Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24-061 

Jtb Yonsei University, Seoul 
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657 million BB pairs collected with the Belle detector at the KEKB asymmetric-energy collider. 
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In the standard model (SM), a mostly left (right-- 
handed photon emitted from a B (B°) meson is ex- 
pected in the b —> s^/ transition. Hence a small time- 
dependent CP asymmetry is predicted in decays of the 
type B — > jcpl where /cp is a CP eigcnstate. New 
Physics (NP) may lead to deviations from the SM ex- 
pectation by introducing different photon polarizations 
in the transition, and can be probed via experimental 
measurements of CP violation parameters Q. In multi- 
body final states B° -> P°Q°j with P° and Q° being C 
eigenstates, the same argument holds 0]. Measurements 
of the time-dependent CP asymmetry in B° — > K^tt ^ 
have been reported by Belle and BaBar based on 535 
and 431 million BB pairs [3, HI, respectively. In this pa- 
per, based on 657 million BB pairs collected with the 
Belle detector Q at the KEKB asymmetric-energy e + e~ 
collider Q, we report the measurement of CP- violation 
parameters on a new channel B° — > Kgp°j where the 
B° decay vertex can be reconstructed from two charged 
pions from the p° decays. 

At the KEKB, the T(4S) is produced with a Lorentz 
boost of f3j = 0.425 along the z axis, which is defined as 
the direction antiparallel to the e + beam direction. In 

the decay chain T(45) — » B°B — » / s i g /tagj where one 
of the B mesons decays at time t S i g to the signal mode 
/si g and the other decays at time t tag to a final state /t ag 

that distinguishes B° and B , the time-dependent decay 
rate is given by: 



P(At,q) 



;-|At|/T B o 

4r R o 



l + q 



S sm(ArrtdAt) 
Acos(Am d At)] (1) 



Here the probability density function (PDF) is normal- 
ized as a function of two variables: B flavor q = +1(— 1) 

when the tagging B meson is B°(B ), and decay time 
difference At between two B mesons. In Eq.(p}, S and 
A are CP-violation parameters, tb» is the B° lifetime, 



Am,d is the mass difference between the two B° mass 
eigenstates. Since the B°B mesons arc approximately 
at rest in the T(4S) center-of-mass system (c.m.s.), At 
can be determined from Az, the displacement in z be- 
tween the / S ig and /tagdecay vertices: At ~ Az/((3"fc). 

The Belle detector [6j is a large-solid-angle magnetic 
spectrometer that consists of a silicon vertex detector 
(SVD), a 50-layer central drift chamber, an array of aero- 
gel threshold Cherenkov counters, and an electromag- 
netic calorimeter (ECL) comprised of CsI(Tl) crystals 
located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux-return located 
outside the coil is instrumented to detect AT° mesons and 
identify muons. 

High energy photons are selected from isolated ECL 
clusters with no corresponding charged track and c.m.s. 
energy satisfying 1.4 GeV < E* < 3.4 GcV. The photons 
are also required to lie in the barrel region of the ECL, 
and have a photon shower-like shape E9/E25 > 0.95, 
where Eg and E25 are the energies summed in 3 x 3 and 
5x5 arrays of crystals around the center of the shower, 
respectively. To reduce background from ir° — > 77 or 77 — > 
77, a 7T° (77) veto is applied with C^o < 0.25 (£,, < 0.2), 
where C„ou,) is a 7r°(?7) likelihood described in Ref. 0. 

Neutral kaons {Kg) are reconstructed from two oppo- 
sitely charged pions whose invariant mass lies within 15 
MeV/c 2 of the Kg nominal mass. Requirements on im- 
pact parameter and vertex displacement are applied || . 

Charged tracks are required to originate from the vicin- 
ity of the interaction point (IP). Charged pions should 
have kaon and pion identification likelihoods consistent 
with the pion hypothesis. This requirement has an effi- 
ciency of 85% with a 7% kaon fake rate. 

The decay B + — > A' + 7r~7r + 7 is also reconstructed to 
study the A7T7t system and to serve as a control sample. 
Charged kaons are selected from charged tracks and re- 
quired to be identified as a kaon with 86% efficiency and 
an 8% pion fake rate. 

We form two kinematic variables: the energy difference 
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AE 



E* 



E* 



and the beam-energy constrained 

is the 



) 2 where Et 



maSS M bc = V^bcam) 2 - (PB) > "bean, 

beam energy in the c.m.s., and E* B and p* B are the c.m.s. 
energy and momentum of the reconstructed B candi- 
date, respectively. The requirement — O.lGeV < AE < 



0.08 GeV is applied. The 



r + and K 



Stt 7r + invari- 



ant masses are required to be less than 1.8 GeV/c 2 . 

The B° — > Kgp°-f candidates are selected from the 
Kgir + ir~j sample by requiring the ir + ir~ invariant mass 
to lie in the p° region, 0.6GeV/c 2 < m nv < 0.9GeV/c 2 . 
Since the p° is wide, other modes, such as A'* + 7r~7 
may also contribute. We first measure the effective CP- 
violating parameters, S e g and A e g, using the final sample 
and then convert them to the CP- violating parameters 
of B° — > Kgp°~f using a dilution factor T>, which is dis- 
cussed later. 

In order to suppress the background from light quark 
pair production qq (e + e~ — > qq with q = u,d,s,c), 
the selection based on an event likelihood ratio 7Z = 
£ s i g /(£sig + jCbkg) is applied. The likelihood for sig- 
nal (£ S ig) and background (£bkg) is formed by combin- 
ing a Fisher discriminant T that uses extended modi- 
fied Fox- Wolfram moments [l(| and the polar angle of 
the B meson in the c.m.s. (cos#b). If there are multi- 
ple candidates in an event, we choose the candidate that 
has the largest 1Z. We define AE sideband events with 
the criteria 0.1 GeV < AE < 0.5 GeV and by vetoing 
B — > (Ktt )j. Other backgrounds that pass our selec- 
tion criteria are classified as (a) B — > K*j background, 
(b) B X s -f background, (c) BB background, which in- 
cludes a generic b — > c component and charmless P-decay 
component. 

The self-cross-feed (SCF) component consists of events 
in which not all tracks are from the signal side, in contrast 
to a true signal with all tracks correctly assigned. We 
define the total signal yield as the sum of SCF and true 
signal events. The fraction of SCF in the total signal 
yield ranges from 5.5% to 10.8% in the Mb c signal region 
defined as 5.27 GcV/c 2 < M hc < 5.29 GcV/c 2 , depending 
on the kaonic resonance. 

We obtain the flavor q, tagging quality factor r 6 [0, 1], 
and signal and tag side vertices from the procedure de- 



scribed in Ref. [ill], using the two charged pions to re- 
construct the signal vertex position. There is no flavor 
discrimination when r = 0, while the flavor tagging is 
unambiguous when r = 1. Events are sorted into seven 
r-bins. 

We obtain 299 events in the signal Mhc region after 
vertexing. An unbinncd maximum likelihood fit to Mhc 
is applied with signal and qq yields floated in each r-bin. 
In the signal Mh c region, we find 212±17 total signal with 
a fraction of 6.0% SCF, 53.4 ±2.6 continuum, along with 
7.8 if*7, 21.5 other X s ^, and 9.0 BB events. The true 
signal shape is parameterized as a Crystal Ball function 
with a width corrected using the B + — > A' + 7r~7r + 7 data. 
The qq shape is obtained from the AE sideband events. 




M bc (GeV/c 2 ) 

FIG. 1: Mbc distributions for B° — > K^tt + tv^j events. Points 
with error bars are data. The curves show the results from the 
r dependent M bc fit. The dashed and dash-dotted curves are 
the qq and all BG. The thin curve is the total signal including 
SCF and the thick curve is the total PDF. 



For other backgrounds and SCF, shapes and r-dependent 
fractions are obtained from MC since their contributions 
are limited. Figure [T] shows the Mh c projection of the fit. 

The ScS and A c s parameters are extracted from an 
unbinncd maximum likelihood fit to the At distribution. 
The likelihood function for each event is: 



/-t-oo 
d(At') 
-OO 

+ f ol P oi (Ati), 



o 

(2) 



where j runs over the signal, SCF and four BG compo- 
nents (qq, -fsT*7, other X s ^, BB). The fraction of each 
component (fj) is calculated using the r dependent Mb c 
fit result on an event-by-event basis. Rj is the resolution 
function in the At. 

The PDF for the signal At distribution P s ig, is given by 
a modified form of Eq.([T]), which incorporates the effect of 
incorrect flavor assignment. The parameterization of R 
is the same as the one used in the B° — > <f> K° analysis 0] . 
The PDF for SCF (Pscf) is the same as for the signal 
except for using a shorter lifetime 1.16 ± 0.02 ps, deter- 
mined from the MC study. The same functional forms 
for the PDF and resolution are used for the K*^f, other 
V s 7 and BB components but with other lifetime val- 
ues obtained from MC and CP parameters fixed to zero. 
The PDF for qq background events, P q q, is modeled as a 
S function convolved with a double-Gaussian resolution 
function R q q. The parameters in R q q is determined from 
a fit to AE sideband events. P Q \ is a Gaussian function 
that represents a small outlier component with a fraction 

/ol- 

The only free parameters in the CP fit to B° — ► 
Kgir + ir~j are iS c ff and A c s, which are determined by 
maximizing the likelihood function L = FT i P; for events 
in the Mh c signal region. We obtain S e g = 0.09 ± 
0.27(stat.)lo.o7(syst.) and A cS = 0.05 ± 0.18(stat.) ± 
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FIG. 2: Fit projections on the At distributions with (a) 
q = +1 and (b) q — —1 for events with r > 0.5. The solid 
curves are the fit while the dashed curves show the back- 
ground contributions. The raw asymmetry as a function of 
At is shown in (c) with a fit curve superimposed. 



0.06(syst.). We define the raw asymmetry in each At 
bin by (N+ -N-)/(N+ + N-), where N+ (iV_) is the 
number of candidate events with q = +1( — 1). Figure [5] 
shows the At distributions and the raw asymmetry for 
events with 0.5 < r < 1.0. 

Various validity checks for our fitting procedure are 
performed. Lifetime fit results for the B + — ► A+7r~7r+7 
and B — > A'g7r+7r~7 modes arc consistent with the nom- 
inal B + and B° lifetimes. No CP asymmetry is seen for 
the control sample B + — > A+7r _ 7r+7. 

We evaluate systematic uncertainties from the follow- 
ing sources. The largest contribution is due to the ver- 
tex reconstruction, where the selection criteria are varied 
to calculate the systematics as A£> c ff =io'oe>i ^A e g = 
±0.03. The <S values for b — > s CP modes arc varied 
from to sin 2<pi to estimate the effect of CP asymmetry 
in the BB background. We fit the data with each fixed 
parameter shifted by its error to evaluate the uncertain- 
ties due to: background and SCF fractions, shapes, At 
PDFs, resolution function, flavor tagging and physics pa- 
rameters Tgo, Amd- Effects of tag side interference are 
evaluated in the same way as in Ref. [3| ■ 

The parameter iS c ff is related to S for Kgp 17 with a di- 



lution factor V that depends on the A"* ± 7r =F components: 



V = 



Sett 

J[\Fa\ 2 



2$t(FXF B )+F* (A) F B (A)] 



f[\F A \* + 2$t(FXF B ) + \F B \ 



(3) 



where Fa , F B are photon-helicity averaged amplitudes 
for B° -> A'g / 9 (7r + 7r-)7 and_B° -> A* ± (A^7r ± )7r = F 7 , 
respectively. The factors F B (K), F B (K) distinguish be- 
tween A* - 7r + 7 and A*+7r~7. The phase space integral 
is over the p° region. 

The charged mode B + — » K + ir~w + 'y is first studied 
using a combination of various kaonic resonances with 
spin > 1 to model the A7T7t system. The amplitude for 
a kaonic resonance A res that decays into particle c and 



subrcsonance r, where r later decays into particles a and 
b, can be modeled by the product of two Breit-Wigncrs: 



M. 



abc\r 



BW(K rcs )BW(r)F K (K res )F r f spin , (4) 



where BW is the relativistic Breit-Wigner lineshapc: 
BW{K rcs ) = l/{m 2 Kl . cs - m 2 abc - im Krcs T Krcs ) and 
BW{r) = l/(m 2 r - m 2 ab - im r T ab ). The width r ab = 

^r(q/qo) 2L+1 f^^J i s a function of q and q , the mo- 
menta of particle a in the r rest frame with mass m a \, 
and m r , respectively. Fx(K res ) and F r are the Blatt- 
Weisskopf penetration factors [l2| for resonances A ros 
and r. The spin factor / sp i n for the resonances used in 

Ap , A*7T, 
= sin# for 



the nominal fit is: / sp m = 1 for A'+(1270) 
A+(1400) -> A*tt with a 1+5 wave; / f 



spin 



A*(1680) -> Kp°,K*n with a l~P wave; / spin = sin0 
for A+(1430) -> K*7r,Kp° with a 2+D wave. Here 9 is 
the helicity angle of subrcsonance r. 

Kaonic resonances with different spin-parity waves do 
not interfere if the decay plane orientation variables are 
integrated out. The total rate is an incoherent sum of 
contributions from different spin-parities with a phase 
space factor. 

To determine the A+p°7 component in B + — > 
A'+7r+7r~7, we study events in a A* region defined as 
\m K+7r - - 0.8961| < 0.075 GeV/c 2 , where most of the 
B + — ► A'+7r~7r+7 signal lies. We fit the rriKmr, m„+„-, 
m K+-K- distributions for events in the signal M^ c re- 
gion. In these fits, the yields and shapes for back- 
grounds are obtained from the corresponding Mb c fit re- 
sults and MC. The fraction and shape for SCF are associ- 
ated with each kaonic resonance from MC. The yields for 
the 2+ AJ (1430)7 components are always fixed based on 
the measured branching fraction B(B — » Aj (1430)7) = 
(1.24 ±0.24) x 10- 5 Jig. 

A two-dimensional (2D) 771^,777^ fit is performed in 
the A* region. Here the SCF is further categorized de- 
pending on whether the A+,7r~ tracks are correctly re- 
constructed. The floating parameters are the rate and 
phase for Ai(1270) -> K*ir relative to Ai(1270) -> Kp. 
In the fit, the yields of 1+, l - kaonic resonances are ob- 
tained from a fit to the tukw-k distribution in the A* 
region. The rates for A+p° relative to A*°7r+ in the 
A* (1680) and A2*(1430) resonances are fixed according 
to the PDG values [lj|, taking into account phase-space 
and isospin factors. The phases for A res — > A*°7r+ rel- 
ative to A rcs -> K+p° in the A* (1680), A 2 *(1430) res- 
onances are fixed to be —0.36 rad [3l and —30° [lij ]. 
respectively. The PDF for the mx-mr distribution in the 
p° region is then obtained. 

The procedure is repeated for different fixed values of 
the phase <j> for Ai(1400) -> A*tt relative to A x (1270) -> 
Kp. For the nominal fit, we use <j> = 15° which gives the 
smallest total x 2 f° r the 2D ttik-k > m-im distribution in 
the A* region and vtlkix-k distribution in the p° region. 
Figure [3fa) shows the distributions for m nn in the whole 
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FIG. 3 



distributions for (a) B + — > A + 7r~7r + 7 and (b) 
B u — > A"s7r + 7r _ 7 . The curves follow the convention in Fig.[T] 
The thin dashed curve is the correctly reconstructed B — + 
Ki (1270)7 signal. 



TABLE I: The number of events for each kaonic resonance 

and various interference terms for the final state As7r + 7r~7. 

'0 „o 



Interf. denotes the interference between K$p' 
= 2&{FXF B ). 



7 and K* + tt 7 





Total 


AV7 


A'* + 7T 7 


Interf. 


F B (K)- 
A s (A) 


A- rea (l+)7 


193.6 


151.0 


35.1 


7.5 


4.4 


(A? (1270)7) 


(167.6) 


(151.0) 


(38.0) 


(-21.4) 


(5.2) 


A rcs (l-) 7 


24.2 


11.3 


8.0 


4.9 


1.3 


A 2 *°(1430) 7 


10.4 


2.2 


6.1 


2.0 


4.5 


Sum 


228.1 


164.4 


49.2 


14.5 


10.2 



events with txikk-k 
m„ n < 0.9GeV/c 2 



ters S Kgp0j = 0.11 ± 0.33(stat.)l^;^(syst.) and Aff 



sidcrcd include: model uncertainty from AT*°(1430)7T7, 
K\ (1400)7, A"i(1270) — > Kui components, uncertainty in 
subresonance rates and phases, SCF uncertainty, isospin 
breaking, the mass and width of kaonic resonances, and 
statistical uncertainty from data. Adding all the uncer- 
tainties in quadrature, we obtain V = 0.83±g J|. 

By combining S c g and the dilution factor T>, we obtain 
SkspO-j = 0-11 ± 0.33(stat.)t°;^( S yst.). 

In summary, we have measured the time-dependent 
CP asymmetry in the decay B° — > Kgp 17 using 
< 1.8 GcV/c 2 and 0.6GcV/c 2 < 
We obtain CP-violation parame- 

0.05/ 
0.09 1 

0.05 ± 0.18(stat.) ± 0.06(syst.). With the present statis- 
tics, the result is consistent with zero and comparable in 
precision to the measurement in B° — > KgTr°j. This is 
the first measurement of CP asymmetry parameters in 
the Ksp°j mode and constrains right-handed currents 
from NP. 
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(China); DST (India); MOEHRD, KOSEF and KRF 
(Korea); KBN (Poland); MES and RFAAE (Russia); 
ARRS (Slovenia); SNSF (Switzerland); NSC and MOE 
(Taiwan); and DOE (USA). 



region. The data distributions for m/t^Tr, mit^, m^w in 
other regions also agree with our model. 

Using isospin symmetry, we assume that the fraction 
and phase of each kaonic resonance channel in the B° 
decay is the same as that in the B + decay. Therefore, 
using the fit results for the kaonic resonance composi- 
tion and substituting the neutral kaonic resonances and 
Kg mass instead of K + mass, we obtain the contribu- 
tions of the terms Aa| 2 , \Fb\ 2 in Eq.Q, as listed in 
Table fl] Figure [3jb) shows the distributions for 

B° -> K a s ir+n-~f. 

Various systematic uncertainties in the dilution factor 
have been investigated. The dominant one comes from 
the modeling systematics. We put additional resonances 
to the default model and repeated the fit procedure. For 
B — > Kaj in which the a has a C parity opposite to that 
of the p°, the dilution factor is: 



V 



J[f_(ao*F-(A > )-f;f + ] 



I\F- 



A+| 2 



(5) 



Here F+ is the photon-helicity averaged amplitude for 
ATi(1270)(-> Allah, while A_ = F A +F B [Eq.©]. In the 
absence of F+, Eq.© reduces to Eq.([3]). The T> value in 
this case is calculated to be 1.00. Other systematics con- 
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